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ABSTRACT 

Previous  application  of  the  Bray-Moss  model  for  premixed 


turbulent  combustion  to  planar,  oblique  and  normal  flames  is  extended 
to  include  the  effect  of  large  but  finite  values  of  the  two  dominant 
characteristic  numbers,  a turbulence  Reynolds  number  providing  a 
measure  of  the  intensity  of  the  turbulence  and  a Damkohler  number 
relating  a turbulence  time  to  a chemical  time.  A classical  perturbation 
analysis  involving  two  small  parameters  proportional  to  the  inverse 
of  these  two  numbers  is  carried  out  to  account  separately  for  the  effects 
of  molecular  transport  and  of  altered  scalar  dissipation  and  for  the 
effects  of  finite  chemical  reaction  rates.  Two  limiting  cases  corres- 
ponding to  highly  oblique  confined  flames  and  normal  or  unconfined 
oblique  flames  are  treated.  Of  particular  interest  in  the  former  case 
is  the  effect  of  the  perturbations  of  the  predicted  orientation  of  the 
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turbulent  reaction  zone.  For  the  unconflned  flames  attention  focuses 
on  the  effect  of  the  perturbations  on  the  turbulent  flame  speed  and  on 
the  change  in  turbulent  kinetic  energy  across  the  reaction  zone.  The 
characteristics  of  the  related  laminr.r  flame  are  introduced  so  that  the 
theory  can  conform  to  the  practice  of  experimentalists  in  correlating 
their  results  for  turbulent  flame  behavior  in  terms  of  such  laminar 
flame  characteristics.  With  respect  to  unconfined  flames,  for  which 
considerable  but  often  contradictory  experimental  data  are  available, 
the  perturbation  analysis  appears  to  yield  qualitative  agreement  with 
a recent  correlation  of  experimental  data  showing  the  effect  of  turbulence 
Reynolds  number.  However,  a comparison  with  results  of  individual 
experiments  is  inconclusive. 


1 


I 

} 


1 

i 

I 


] 

i 


= 

I 

i 

1 

1.  INTl^ODUCTION 

Numerous  experimental  8tudL<!S  of  the  propagation  of  turbulent 
flames  into  premixed  combustible  gas  mixtures  employ  the  properties  \ 

of  laminar  flames  to  correlate  the  observed  turbulent  flame 

speeds.  For  example,  Abdel-Gayed  and  Bradley  [1]  recently  J 

presented  a correlation  of  a large  amount  of  experimental  data, 
involving  many  different  combustible  mixtures,  in  which  the  ratio 
of  turbulent  to  laminar  flame  speed  is  presented  as  a function  of  the 
ratio  of  laminar  flame  speed  to  turbulence  intensity  and  of  a turbulence 
Reynolds  number.  Other  experimental  data,  for  example,  that  of 
Wright  and  Zukoski  [2],  are  insensitive  to  factors  which  are  known 

» 

to  influence  the  laminar  flame  speed.  However,  it  is  beyond  dispute 
that  a regime  exists  in  which  the  turbulent  flame  speed  is  strongly 
influenced  by  the  corresponding  laminar  flame  speed. 

A successful  model  of  turbulent  flames,  having  any  claim  to 
generality,  must  be  able  to  predict  these  effects,  and  must  therefore 
include  the  relevant  fluid  mechanical  and  chemical  phenomena.  On 
the  other  hand,  existing  models  (Borghi  and  Moreau  [3],  Pope  [4] 
and  Bray  and  Moss  [5  , 6]  ) generally  neglect  molecular  transport 
and  effects  of  finite  rates  of  combustion  reactions,  and  treat  only  the 
I limiting  case  where  the  time-average  rate  of  heat  release  is  controlled 


by  turbulent  mixing.  Because  the  properties  of  a laminar  flame  depend 


4 


upon  a balance  between  finite  rates  of  molecular  transport 
and  chemical  reaction,  it  follows  that  comparison  between 
predictions  from  such  theories  and  experiment  is  at  best  inadequate 
and  incomplete. 

I 

1 

I 

Another  problem,  of  considerable  practical  importaace,  is  the  ! 

1 

prediction  of  the  limits  of  stable  combustion.  These  limits  are  known  ' 

to  be  strongly  influenced  by  chemical  kinetic  factors,  and  may  also  be  j 

dependent  on  Reynolds  number,  providing  a further  incentive  to  improve  ' 

current  combustion  models  to  include  these  effects.  | 

The  objective  of  the  present  work  is  to  extend  the  model  of  premixed  j 

turbulent  combustion,  developed  and  exploited  by  the  present  authors 
(Bray  and  Moss  [5,6],  Bray  and  Libby  [7]  and  Libby  and  Bray  [8]),  in 
order  to  include  effects  of  a finite  reaction  rate  and  of  molecular 
transport,  and  to  compare  predicted  trends  with  those  observed 
experimentally.  The  Bray-Libby  and  Libby-Bray  references  will  be 
cited  repeatedly  and  are  therefore  denoted  I and  11. 

Our  presentation  is  organized  as  follows:  we  first  review  briefly 
the  Bray-Moss  model  and  then  discuss  the  fundamentals  required 
for  the  inclusion  of  finite  but  large  values  of  the  two  characterizing 
parameters,  a turbulence  Reynolds  number  and  a Damkohler  number; 
these  are  taken  to  be  infinite  in  previous  applications  of  the  model 


and  accordingly  lead  to  a first-order  theory.  Sections  4 and  5 describe 


I 


a perturbation  analysis  for  the  two  limiting  cases  of  planar  flames 

we  consider  seminal:  highly  oblique  confined  flames  and  normal  or  ! 

unconfined  oblique  flames.  In  the  numerical  computation  of  these  i 

cases,  values  for  certain  parameters  are  needed.  Here  we  employ 

the  values  suggested  by  a comparison  of  the  predictions  of  the  first  ' 

1 

r 

! 

order  theory  with  experiment.  Section  6 relates  the  application  of 

I 

! 

the  perturbation  analysis  to  current  methods  of  presentation  of  j 

experimental  data  and  incorporates  the  relevant  laminar  flame  theory 
in  order  to  provide  estimates  for  several  of  the  parameters  which  are 
thereby  encountered.  Finally,  conclusions  are  drawn  from  the  analysis. 


2.  BRIEF  DISCUSSION  OF  THE  EXISTING  MODEL 

For  completeness  and  to  provide  the  requisite  background  for 
the  developments  presented  here  we  discuss  the  general  features  of 
the  Bray-Moss  model  of  premixed  turbulent  reactions  and  its  application 
to  plane  flames.  Details  are  in  Bray  and  Moss  [5,  6] 
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and  in  I.  The  model  utilizes  a series  of  assumptions;  some  relate 
to  the  aerothermochemistry  of  the  flow  and  are  standard  in  the 
combustion  literature.  Of  central  importance  in  this  regard  is  the 
result  that  the  instantaneous  thermodynamic  state  of  the  reacting 
mixture  depends  solely  on  a reaction  progress  variable  c , in 
pa  rticular 


T = T(c)  = T (1  + Tc) 
o 


(2.  1) 


P = p(c)  = p (1  + Tc) 
o 


(2.2) 


where  T is  the  temperature,  p is  the  mixture  mass  density  and  the 
subscript  "o"  denotes  conditions  upstream  of  the  reaction  zone  where 
c = 0.  The  quantity  T plays  an  important  role  in  the  description  of 
the  phenomena:  from  Eq.  (2.  1)  we  note  that  T = (T  /T  ) - 1 , with 

GO  O 

T the  temperature  downstream  of  the  reaction  zone  where  c = 1 . 

Values  of  T of  practical  interest  are  from  4 to  9.  The  progress  variable 
c can  be  considered  to  be  the  mass  fraction  of  the  product  of  the  one-step 
reaction  normalized  by  its  value  when  reaction  is  complete;  thus  c 
will  usually  be  loosely  termed  the  product  concentration. 

A consequence  of  the  dependence  of  the  instantaneous  state  of  the 
gas  mixture  on  the  value  of  c is  that  the  rate  of  production  of  product, 
i.e.,  the  source  term  in  the  conservation  equation  of  product,  is  also 
a function  of  c alone,  denoted  w(c),  with  a maximum  value  w 

max 
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The  Bray-Moss  model  is  based  on  considerations  of  the 
probability  density  function  for  the  concentration  c . At  a ^iven 
spatial  location  within  the  reaction  /.one  this  pdf  will  involve 
a partition  among  unburnt  (c  = 0)  , burning  (0  < c < 1)  , 
and  all  burnt  (c  = 1)  mixtures.  Strengths  of  OL  , 7 , and  P are 
attributed  to  each  of  these  respectively.  The  assumption  that  the  rate 
of  chemical  reaction  is  fast  and  thus  that  the  rate  of  heat  release  is 
determined  by  turbulent  mixing  implies  that  T « 1 . Accordingly,  the 
pdf  of  the  concentration  c is  dominated  by  delta  functions  at  c = 0,  1. 

The  combination  of  the  dependence  of  thermochemical  quantities  on 
c alone,  together  with  a model  for  the  probability  density  function  for  c , 
permits  quantities  of  fluid  mechanical  interest  to  be  readily  expressed 
in  terms  of  the  mean  product  concentration,  c , with  corrections  of 
order  7 depending  on  integrals  of  the  pdf  denoted  f(c)  within  the 
interior  range  0 < c < 1 . Because  it  only  appears  within  such  integrals, 
details  of  f(c)  are  inessential.  Similar  considerations  apply  to  the 
function  describing  the  chemical  reaction  w(c)  . 

In  I and  II  Favre-averaging  is  employed  (cf.  Favre  [9]) 
because  of  the  resulting  formal  simplification  of  the  describing  equations. 
Thus,  for  example,  we  have 


Pc  c /p  = c(l  - c)-y(l  + tc)M 


w = w y 1- 
max  3 


c w = yw  (I.  - c I_) 
max  4 3 


(2.3) 

(2.4) 

(2.5) 


\ 


I 
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where  M and  I denote  integrals  of  the  sort  alluded  to,  namely 

lx 


-1 

= ^ c(l-c)(l+Tc)  f(c)dc 


(2.6) 


I,  , = ^ c^w(c)  w f(c)  dc  , k - 0,  1 

k+3  wIq  max 


(2.7) 


and  where  {'' ) and  ( )'  denote  mass  averaged  and  fluctuating 

quantities  respectively. 

In  formalizing  the  assumption  of  fast  chemistry,  it  is  convenient 
to  introduce  a Damkohler  number.  Da,  defined  as  a ratio  of  a turbulence 
time  scale  to  a chemical  kinetic  time  scale,  in  particular 


Da  = I_  w t /p  q 
3 max  o o o 


(2.8) 


where  is  a length  scale  characterizing  the  large  eddies  and  is 

the  Favre-averaged  turbulent  kinetic  energy,  i.  e. , q = ^ pu^u^  . Now 
for  fast  chemistry  Da  » 1 but  y Da  = 0(1)  . 

The  formulation  of  the  conservation  equations  based  on  these 
assumptions  focuses  on  the  equations  for  turbulent  kinetic  energy,  q ; 
the  mean  product  concentration,  c ; and  the  intensity  of  the  product 
fluctuations,  Pc  'c"  . A closure  problem  exists  so  that  some  terms 
must  be  modelled.  At  the  present  time  a gradient  approximation  with 
an  eddy  transport  coefficient  based  on  the  Prandtl-Kolmogoroff  model 


In  the  course  of  the  analysis  the  tilde  is  also  introduced  to  denote 
several  new  variables.  No  confusion  should  result. 
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is  employed  with  a modification  to  account,  at  least  roughly,  for  variable 
density  effects.  We  have  (cf.  [8]) 

= aq^/j(p/p^)P  (2.9) 

where  a is  a constant,  is  a scale  of  the  large  eddies  within  the 

reaction  zone,  and  p is  an  exponent  to  be  assigned.  With  Eq,  (2.  9) 
a typical  flux  term  becomes 


puj  g"/p  = - 


(2.  10) 


In  the  calculations  of  I and  II  molecular  transport  is  neglected, 
implying  consideration  of  an  infinite  turbulence  Reynolds  number. 


R 


T 


= t tv 


(2.  11) 


where  V is  the  molecular  kinematic  viscosity. 

There  remain  to  be  considered  the  dissipation  of  turbulent  kinetic 
energy  and  of  the  concentration  fluctuations.  With  respect  to  the  former, 
denoted  tP  , the  nature  of  the  equations  indicates  that  within  the  thickness 
of  the  reaction  zone  the  dissipation  of  turbulent  kinetic  energy  must  be 
negligible  as  is  found  to  be  the  case  when  turbulence  undergoes  rapid 
distortion  (cf.,  e.g.,  Batchelor  [lO]).  On  the  contrary  the  annihilation 
of  concentration  fluctuations  by  molecular  processes  is  found  to  be 
important  and  is  modeled  close  to  the  usual  fashion,  but  with  variable 
density  effects  included,  at  least  roughly:  we  take 

X = Pc"c"n^ 


(2. 12) 
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where  JL^  is  another  length  scale  attributed  to  scalar  dissipation  and 
where  k is  another  exponent  to  be  assigned. 

The  ecuations  which  result  from,  these  considerations  are  specialized 
and  subsequently  applied  to  planar  turbulent  flames.  In  the  case  of 
normal  flames  and  unconfined  oblique  flames  no  additional  assumptions 
are  required;  in  the  case  of  confined,  oblique  flames  of  specified 
orientation  identified  with  the  angle  6,  6=  90°  for  normal  flames,  an 
additional  approximation,  namely  that  the  mean  streamline  is  undeflected, 
results  in  significant  simplification  and  yet  is  physically  reasonable  for 
confined  turbulent  flames,  e.g.  , as  in  Wright  and  Zukoski  [2].  A con- 
sequence of  this  assumption  is  that  there  results  an  explicit  expression 
for  the  Reynolds  stress,  i.  e.  , pu"v"/p  in  Favre-averaging,  in  terms 
of  the  mean  product  concentration.  An  important  physical  implication 
which  derives  from  this  consideration  is  that  the  acceleration  through  the 
reaction  zone  of  a confined  oblique  flame  is  due  to  the  gradient  of  Reynolds 
stress,  an  implication  not  based  on  a gradient  assumption  but  rather  on 
the  assumed  kinematics  of  the  flow.  This  contrasts  with  a normal 
flame  in  which  the  pressure  gradient,  although  thermodynamically 
negligible,  accelerates  the  gas.  The  difference  in  mechanism  in  the 
two  cases  suggests  that  the  appropriate  models  of  the  laminar  flame 
within  the  turbulent  reaction  zone  used  to  estimate  f(c)  may  likewise 
be  different. 

Although  for  brevity  we  characterize  flames  as  either  oblique 
or  normal,  these  two  cases  are  considered  to  correspond  respectively 


to  those  with  constrained  and  unconstrained  mean  streamlines,  since  an 
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oblique  flame  involving  a constant  tangential  velocity  and  thus  a curved 
mean  streamline  corresponds  exactly  to  a normal  flame  insofar  as  our 
analysis  is  concerned. 

The  orientation  of  the  confined  flame  determines  the  relative 
importance  of  two  competing  effects  related  to  the  interaction  of  the 
turbulence  and  the  heat  release  associated  with  chemical  reaction. 

If  the  flame  is  nearly  normal,  dilatation  dominates  and  the  turbulent 
kinetic  energy  downstream  of  the  reaction  zone  is  less  than  that  in 
the  oncoming  stream.  On  the  contrary,  for  oblique  flames  the  generation 
of  turbulent  kinetic  energy  by  interaction  of  the  Reynolds  stress  and 
the  gradient  of  the  mean  velocity  dominates  and  the  turbulent  kinetic 
downstream  exceeds  that  upstream.  This  situation  leads  to  an  important 
special  case,  termed  strong  interaction,  in  which  the  downstream 
kinetic  energy  overwhelms  that  in  the  oncoming  stream. 

The  resultant  mathematical  formulation  leads  to  a double  eigenvalue 

problem:  the  turbulent  flame  speed  in  the  form  ^ /'a^  and  the  ratio 

o o 

of  two  turbulent  kinetic  energies.  Q = q /q  , are  predicted  for 
normal  flames.  In  the  strong  interaction  case  the  ratio  of  the  turbulent 
flame  speed  to  the  turbulent  kinetic  energy  downstream  of  the  flame, 
orientation  of  the  flame  B are  given  as  part  of  the 
solution.  Thus  the  mathematical  formulation  corresponds  to  a more 
complex  version  of  that  arising  in  laminar  flame  theory  in  which  a 
single  eigenvalue,  involving  chemical  kinetic  and  fluid  mechanical 
quantities,  is  predicted. 


In  I and  II  considerable  attention  is  devoted  to  comparison 
between  prediction  and  appropriate  experimental  results.  As  indicated 
earlier  this  comparison  is  limited  for  a variety  of  reasons  including 
the  limitations  of  the  theory.  Nevertheless,  in  II  a comparison  between 
the  orientation  of  flames  in  the  strong  interaction  limit  suggests  that 
the  two  exponents  incorporating  at  least  roughly  variable  density  effects 
in  the  modelling  in  the  form  (p  + k)  should  equal  two;  the  implication 
from  this  result  is  that  the  usual  modelling  should  be  altered  when 
significant  density  variations  occur. 


3.  INCLUSION  OF  FINITE  DAMKOHLER  AND  REYNOLDS  NUMBERS 
Wo  now  consider  the  extension  of  the  previous  application  of  the 
Bray-Moss  model  to  planar  turbulent  flames  so  as  to  include  finite 
values  for  the  two  characteristic  parameters,  the  Damkohler  and 
Reynolds  numbers.  The  flow  we  consider  is  shown  schematically  in 
Fig.  1.  The  Fa vre-averaged  balance  equations  incorporating  the 
assumptions  discussed  earlier  are  as  follows: 


Continuity; 


p u = 0 u 
o o 


Species: 


(3.  1) 


— dc 


w 


(3.2) 
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Species  fluctuation: 


N It  He 


Pu  — ) -2PC  u ;t-  » 2c 


(pc  c u ) + 2—  (pDc  ) - X 


(3.3) 


Turbulence  kinetic  energy: 


dq  //  ff  du  _ ,/  //  dv 

P a = - p\i  \i  — -puv  — 
dx  dx  dx 


d ,1  _ « » » V . d , //  , . — 

- — (?  P u V V ) + — (v  f , ) - V 
dx  a Of  dx  ct  Ofl 


(3.4) 


To  complete  these  equations  the  modelling  of  the  turbulence 
terms  alluded  to  earlier  must  be  incorporated.  In  addition,  for  present 
considerations  the  terms  associated  with  the  molecular  transport 
coefficient,  pD,  and  with  the  viscous  stress  tensor 


d V dv. 


Z.  2.  * , / 0?  , 0 ' 

3 dx-  Of0  ^'dx„  dx  / 


(3.5) 


and  the  effect  of  finite  Reynolds  numbers  on  the  scalar  dissipation  must 

be  incorporated.  In  doing  so  and  in  anticipation  of  a perturbation  treat- 

- 1 1 

ment  we  introduce  a small  parameter  €=  {a.  a t /v  ) =(aR  )"*• 

o o o X ’ ’ 

a second  perturbation  parameter  will  be 


^ -1 

(5  = (a  I w i /p  q ® ) 

3 max  o o ^ o 


(3.  (.) 
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To  treat  the  extra  terms  we  make  several  additional  assumptions: 

i)  The  molecular  Schmidt  number  is  taken  as  unity  so  that  pD  = p . 

ii)  The  product  pp  is  taken  as  constant  so  that  pp  = > 

2 ~ ' ~2  2'-~\ 
y = P (1  + Tc)  , and  P = P , (1  + Tc)  + T c (1  - c)  ) + 0(y)  . 
o o . / 

iii)  The  scalar  dissipation  is  modified  to  include  a Reynolds  number 
correction  so  that  (cf.  Eq.  (2.  12)  ) 


,T/.  >k 


X =C{p/p/q^  Pc'c"'^!  +(CjP/aq^£j)/£2 


(3.7) 


iv)  The  correlation  c*(dc*/3x)  is  replaced  by  id/dx{pc"c  /p)  with 
the  consequence  that 


dc^  — ~dc  d H 1 »d/Pcc\ 

pD^^  pP  — +T-pyc  -tp  V 7—  f ^ - ■ ) 

dx  dx  dx  o o dx  \ ~ / 


But  pPc"  = p P c^ll  - c)  T(T +2) +0(y) , pc'c* /p  = c(l  - c)  + 0(y) 
o 

so  that 


q t 

Sc  ^ T ..  - O O ^ 


pD^-pp^c^;3jy;(l  + rc)P  ^ 


+ t(m1  -2c)-^i^^^)  (T  +2)  - id  -2c)  (1  +TC)]  j ^ 
' " 1 +TC  ' / ^ 


— _ _ dc 

Pl-^.RE,  - 


(3.8) 


where  R = (q^  and  where  E,  = E,(c)  is  implicitly 

o o I 11 

defined.  Two  contributions  to  this  mean  molecular  transport 
term  can  be  identified;  the  first  is  associated  with  the  mean 


IS 


viscosity  coefficiRnt,  the  second  with  fluctuations  in  that  coefficient. 
Both  can  be  readily  incorporated  with  the  single  reasonable 
assumption  cited  here. 

v)  The  correlation  c*c*(Sc'7ax)  is  replaced  with  ( i/3)d/dx  (pc"' ^/p ) 
with  the  consequence  that 


pDc"  ~ = pVc" 
ox 


dc 

dx 


ti  — Pl/C 


ff  2 


But  pVc'^^pV  c (1 -c  ) fl +2T  + 
o \ 

c(l-c)(l-2c)  + 0(7)  so  that 


Jl  i'P£ll 

dx  ^ — 

P 


) 


. t(t  + 2)c  ) + 0(y),  pc"^/p  5= 


/»  Sc  ^ 
pDc  — - 
ox 


pp^eR(i+rc)' 


/-w/  / 5 

c ( 1 - c)  T 


(r  + 2) 


i 1+2T+T  -t(t+2)c  , 1 ,,  ~ 

* ^ ■ + T ( 1 +TC 

(1+Tc)  ^ 


/ dx 


(3.9) 


where  ^2^^^  implicitly  defined.  Note  again  that  we  are  able  to 

incorporate  without  difficulty  the  molecular  transport  of  fluctuations 

in  concentration  due  to  the  mean  and  fluctuating  viscosity  coefficient. 

Vi)  The  description  of  the  flux  of  turbulent  kinetic  energy  due  to 

molecular  processes  as  contained  in  the  term  v *f 

<X  al 


is  more 
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difficult  to  rationalize.  However,  we  incorporate  both  the  effect 
of  fluctuations  in  the  viscosityr  coefficient  and  the  transport  of 
turbulent  kinetic  energy  by  the  mean  molecular  coefficient. 
Without  full  justification  we  take  at  least  provisionally 


du 


V f •,  = 7 pVa  ^ + pvv 

a Ofl  3 dx 


dv  . d ~ 


(3.10) 


The  effective  transport  coefficients  pUa"  and  puv"  can  be 
computed  by  application  of  the  approach  in  1 1 involving  self- 
consistent  estimates  for  the  velocity  components  of  packets 
of  unburned  and  fully  burned  gases  within  the  reaction  zone. 
We  find 


ppu"  = - p p y T(r  + 2)  ^ + 0(y) 
To  dx 


pP  P T(T  + 2)  ,~ 

pPv*  = — , ^ + OO") 

tan  9 dx 


so  that 


a al 


- PV^€R 


2 

tan  9 


UoPtT^T  + 


~ .2 


2) 


Vdx  / 


JL  ^ \ 

Po  dx/ 


(1  + TC)*" 


(3.11) 


In  order  to  treat  the  flux  of  turbulent  kinetic  energy  as  thoroughly 
as  that  of  mean  concentration  and  mean  intensity  we  would  need 
an  extension  of  the  Bray-Moss  model  to  a joint  velocity-concentration 
probability  density  function. 
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vLL)  The  viscous  dissipation  is  neglected  even  at  finite  Reynolds 
numbers. 

With  these  preliminaries  we  are  prepared  to  proceed  with  the 
treatment  of  the  effects  of  finite  Damkohler  and  Reynolds  numbers. 


4.  APPLICATION  TO  THE  CASE  OF  STRONG  INTERACTION 

We  apply  the  equations  which  result  from  imposition  of  the 
several  considerations  discussed  earlier  to  the  case  of  oblique,  planar 
turbulent  flames  under  conditions  such  that  the  turbulent  kinetic  energy 
generated  by  shear  within  the  reaction  zone  overwhelms  that  in  the  approach- 
ing stream.  In  this  case,  which  is  termed  strong  interaction,  it  is  convenient 
to  replace  q by  Q = (Q  - Q )/( 1 - Q ) , where  Q = q/q  . The  two 

eigenvalues  which  are  to  be  determined  as  part  of  the  solution  and  which 
provide  predictions  of  physical  interest  are 


IC=  u^/q 
o ' 


2 

tan  0 


(4.2) 


where  c^  = I^/l^  . These  two  eigenvalues  determine  the  turbulent  flame 
speed  as  a multiple  of  the  turbulent  kinetic  energy  downstream  of  the 
flame  and  the  orientation  of  the  flame  with  respect  to  the  oncoming  flow. 
For  this  case  Eq.  (3.  2)  becomes 


_1  ( 
dx  V 


P u 
o o 


yw  L 
max  3 


P u 
o o 


(4.3) 


3 


\ 

i 

f 
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which  will  prove  useful  for  subsidiary  calculations.  If  £q.  (4.3)  is 

used  to  eliminate  the  term  involving  w from  Eq.  (3.  3), 

max 

we  have  the  equation 


PK 


T dP^  dP  _ 

dx  ^ ~ dx  dx  ^m 

P “ 

o o 


dc^ 


/ ^ _ _±  ( ^ T 

\dx  dx\^~  dx.  ^ 
' p u 


o o 


Pu  (1  -Q)c  (1  -c) 


:i-- 


p ' 


, , ■ 1 +m  ~ . 

P^(l  +Tc)  c(l  - c ) 


O'  ' A II 


P u 
o o 


^ d / dc  '\ 

A II  ' 


P u 
o o 


(4.4) 


Equation  (3.4)  becomes 


-1,'.^  dQ'  /C(l  -Q)T  N ^ 

dx  ~ dx  / dx  ^ 1 . ~ T c 

p u 1 + Tc 


>'  dx 


o o 


dx  . 

P u 

o o 


r((|  tan^e  + l)— ^ 


KP  T (t  + 2) 


Sf-f#)..-) 


(4.5) 


where  P = yfl  +tc)M; 

m = p + k : and 
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dG  _ G ^ . KT^c)  - C 4 ( «(  (l  + " * 

dc  ® ^ 1 +TC  dc  \ ' 


T^(T  +2)  (1  +Tc)  + 


Xo),..re,-) 


d Q _ G 
dc 


(4.9) 


(4. 10) 


The  perturbatioa  equations 

Equations  (4.7)  -(4.  10)  are  in  a form  suitable  for  a perturbation 
analysis.  The  appropriate  expansions  are  as  follows: 


Q=Q  +€Q,  + 6Q,  + ... 

o 1 2 

S=S  +€S,  +6S,  + ... 

o 1 2 

G = G +eG,  +6G,  + ... 

o 1 2 

P = 6(P  + . . . ) 

o 

D = 6(D  + ...  ) 

o 

p = 3 + C3,  + 6P.,  + . . . (4.11) 

o 1 Z 

*■  = K +fK,+6K,  + ... 
o 1 2 

y = 6(y  + ...  ) 

o 


c ~ c 4 . . . 
m mo 

M = M + . . . 

o 

R - R + . . . 
o 

Although  not  required  al  this  staK**  in  the  analysis,  the  x-coordinate 
must  also  be  considered  expanded  accordinji  to  x(c)  = ’‘,j(^)  "*■  (C  ) 4 
Sx^fc)  4 . . . . 
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Note  that  not  all  of  these  expansions  are  independent  but  they  are 
formally  consistent  as  shown. 

When  Eqs.  (4.  11)  are  substituted  into  Eqs.  (4.7)-(4.  10), 
the  terms  independent  of  ( and  6 yield: 


which  are  the  equations  for  the  strong  interaction  case  treated  in  I and 
II  when  c is  taken  as  the  independent  variable.  The  solutions  are 
subject  to  boundary  conditions  which  we  shall  u.scuss  subsequently 

A 

and  yield  as  eigenvalues  P and  K . 

o o 

Collection  of  the  first  order  terms  in  e and  6 in  Eqs.  (4.7)  - 
(4.  10)  yields  the  following  equations,  in  which  for  compactness  we 


( 


I 


J 


employ  the  Kronecker  delta  , with  j = I corresponding  to  the 
C -perturbation  and  j = 2 to  the  6 -perturbation: 
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dP 
c 

dc 

dG. 

1 

dc 


(4.  16) 


G , G.  S.  , /€  U.  T _ 

(T-r)=  77r-o^^~(7^) 


+ 1 u • 

/ o 


t‘'(T+2)(1  +Tc)S^  + ^ G )(1+Tc)^| 

o / 


dQ.  G , G.  S.  ^ 
I _ _0  _l_  \ 

~ S V G " S / 
dc  o o o 


(4.  17) 


(4.  18) 


For  the  perturbation  solutions  corresponding  to  finite  Reynolds  numbers, 
i.e.,  to  e,  i = 1 , Eqs.  (4.  15),  (4.  17),  and  (4.18)  are  three  equations 
for  Sj  , Gj  , and  . For  the  perturbation  solutions  corresponding 
to  finite  Damkohler  numbers,  i.e.,  to  6,  i = 2,  we  must  evaluate 
and  . Equation  (4.  3)  and  the  definition  of  P lead  to 


q 

/ o ^ 


(1  + Tc)  q_ 


(4. 19) 


where  ^*  = (aCA,/(2c  “ ^ ^ 

\ 1 mo  2 / 


taken  to  be  constant.  The  parameter 
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appears  in  I and  11  and  contains  several  empirical  constants 

evaluated  from  the  phenomenology  of  turbulence  for  constant  density 

flows.  Equations  (4.  19)  and  Eq.  (4.  16)  determine  D so  that  the 

o 

left  side  of  Eq.  (4.  15)  can  be  expressed  in  terms  of  the  first  order 

solution.  Thus  again  Eqs.  (4.  15),  (4.  17)  and  (4.  18)  become  three 

equations,  in  this  case,  for  , and  Q^. 

Note  that  the  last  factor  in  Eq.  (4.  19)  involves  explicitly  the 

ratio  of  the  turbulent  kinetic  energies  downstream  and  upstream  of 

the  reaction  zone,  a ratio  large  compared  to  unity  in  the  strong  interaction 

case  being  considered  here.  However,  we  avoid  specifying  this  ratio 

a priori  by  introducing  a new  parameter,  1?l  = q /q  ) . 

o o ® o 


The  boundary  conditions 

In  terms  of  the  spatial  coordinate  x,  boundary  conditions  are 

imposed  at  which  correspond  in  terms  of  ^ to  0,  1 . In  addition 

our  basic  formulation  involves  the  well-known  "cold  boundary  problem" 

which  requires  for  its  resolution  specification  of  a value  of  c = c ; 

o 

for  c < c^  no  chemical  reaction  occurs.  ' In  the  present  context  this 
condition  corresponds  to  specification  of  a mean  ignition  temperature. 

In  I it  is  shown  that  for  the  highly  oblique  flames  associated 

with  the  strong  interaction  case,  the  characteristics  of  the 

flame  are  essentially  independent  of  c^  , a physically  satisfying  result. 

sj; 

This  condition  implies  that  at  c = c^  the  dependent  variables 
are  continuous  but  with  "jumps"  in  their  derivatives. 


i 
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Since  the  describing  equations  both  for  the  first  order  and 

perturbed  solutions  can  be  readily  solved  for  0 s 5 « 1 our 

o 

boundary  conditions  effectively  apply  at  ? , 1 . 

o 

Consider  the  solutions  in  the  range  0 ^ ^ ^ ^ « 1 ; if  the 

o 

right  side  of  Eq.  (4.3)  is  set  to  zero,  we  find 


S = c 
o 

and  subsequently  from  Eq.  (4.  15) 


(4.20) 


(4.21) 

(4. 22) 


Furthermore,  Eqs.  (4.  13)  and  (4.  14),  simplified  as  appropriate  for 

c ^ 0 , Q 1 , yield 
o 

(4.23) 


Similarly,  if  Eqs.  (4.  17)  and  (4.  18)  are  combined,  quadrature  leads 
to  the  result 


G. 

1 


i = 1 


(4.24) 


= 0,  i = 2 


(4.25) 


provided  it  is  recognized  that  R ~ 1 . 

o 

Equations  (4.20)  and  (4.  23),applied  at  c = c^,  provide  the  boundary 
conditions  for  the  first  order  solutions  already  presented  in  I and  II  and 


determine  the  ei>;onvalue8  (5^^  and  • Likewise,  as  we  shall  see, 
the  pairs  of  equations,  r'(|s.  (4.21)  and  (4.24)  and  (4.22)  and  (4.25), 
determine  the  perturbation  eigenvalues. 

The  boundary  conditions  to  be  imposed  at  the  other  end  of  the 

range  of  c , i.  e.,  in  the  neighborhood  of  c = 1 , are  more  complex; 

their  detailed  consideration  is  discussed  in  Appendix  1.  For  present 

purposes  it  is  sufficient  to  note  that  the  asymptotic  behavior  of  the  first 

order  and  perturbation  solutions  as  c ~ 1 is  such  that  no  arbitrariness 

is  involved.  Accordingly,  a numerical  integration  initiated  in  the 

neighborhood  of  c = 1 can  be  carried  out  for  decreasing  c with 

the  boundary  conditions  at  c = c^  satisfied  by  appropriate  selection 

of  the  eigenvalues.  The  numerical  analysis  of  the  first  order  solutions 

requires  iteration  since  the  asymptotic  solutions  as  1 

involve  the  eigenvalues.  However,  the  perturbation  solutions 

can  be  generated  in  a form  such  that  only  a single  integration 

from  the  neighborhood  of  c ~ 1 to  c^  ~ ^ is  required. 

o 

Treatment  of  the  perturbation  equations  Involves  specification 
of  a form  for  the  quantity  which  is  the  ratio  of  products  of  a 
length  scale  and  the  square  root  of  the  turbulent  kinetic  energy.  This 
situation  is  in  contrast  to  that  prevailing  for  the  first  order  solutions; 
in  terms  of  the  mean  product  concentration  these  can  be  determined 
without  specification  of  a length  scale  model.  Here  for  simplicity 
and  with  the  expectation  that  our  conclusions  will  not  be  affected  by 
more  complex  models, we  take  1 which  implies  that  the  length 
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scale  is  inversely  proportional  to  the  square  root  of  the  turbulent 
kinetic  energy. 


The  form  of  the  perturbation  solutions 

In  addition  to  the  parameters  which  determine  the  first  order 

solutions,  namely  t , C,  m,  and  c^  , we  must  specify  for  the 

perturbation  solutions  the  individual  contributions  to  m,  i.e.,  p and 

k,  the  parameters  ♦ and  c , and  the  coefficients  C,  and  Tn  . 

mo  1 o 

However,  it  is  possible  to  carry  out  the  numerical  analysis  of  the 

perturbation  solutions  so  that  for  given  first  order  solutions  only  one 

integration  is  required  without  a priori  specification  of  Cj  and  7?!^. 

To  illustrate  consider  the  solutions  for  Q.  in  the  form 

1 


Q.=6.,Q 

I lI  ip 


+ 


6._  Q.,  , i 1,  2 (4.26) 

i2  o i3 

where  Q.  --  Q.  (c),  n = p,  1,  2,  3 provide  the  constituent  parts  of 
in  in 

the  perturbation  solution  Q^(c).  Similar  expressions  prevail  for 
and  G.  . When  these  forms  are  substituted  into  Eqs.  (4.21)  and  (4.24) 

i ^ 

for  i = 1 and  into  Eqs.  (4.22)  and  (4.25)  for  i = 2 with  c = c^  , we 

obtain  for  each  perturbation  two  algebraic  equations  for  the  two 

pairs  of  eigenvalues  in  the  form  (P./P  ),  ) . Values  for  C, 

1 o i o 1 

and  are  required  if  specific  values  for  the  perturbation 
eigenvalues  are  desired.  However,  the  solutions  to  the  algebraic 
equations  can  be  arranged  to  yii'ld 


Z7 


'>,2 


?,/P  = »!  b,, 
2 o o 21 


K_/k 

2 o o 22 


(4.27) 

(4.28) 

(4.29) 

(4.30) 


where  the  and  b_  coefficients  are  known  and  depend  only  on 

the  first  order  solution,  on  the  parameters  p and  k , and  on  ^ 

‘^mo  ‘ possible  to  assess  readily  the  sensitivity  of 

the  perturbation  eigenvalues  to  the  extra  parameters  C and  JW 

1 ' o 

The  form  of  Eqs.  (4.  '7)  - (4.30)  warrants  comment.  The 
P'-rcurbation  eigenvalues  corresponding  to  large  but  finite  values 
of  the  Reynolds  number  are  due  to  two  separate  effects:  the  a^^ 
and  aj2  coefficients  arise  from  the  molecular  transport  terms  in 
the  conservation  equations.  A second  effect  is  contained  in  the  b^^ 
and  bj2  coefficients  and  is  associated  with  the  parameter  and  thus 
with  the  influence  of  Reynolds  numbers  on  the  scalar  dissipation. 

The  situation  regarding  the  influence  of  finite  but  large  Damkohler 
numbers  is  different  and  is  seen  to  be  associated  with  the  parameter 
alone.  According  to  the  present  analysis  7n^  is  proportional  to 
the  ratio  of  turbulent  kinetic  energies,  (q„/q^)  with  the  proportionality 
factor  depending  on  , i.  e. , on  an  integral  of  the  distribution  of 
product  through  the  reacting  surfaces  in  the  form  f(c). 
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The  physical  significance  of  the  perturbation  eigenvalues 

Of  primary  interest  from  the  numerical  results  is  the  effect  of 
large  but  finite  values  of  the  Reynolds  and  Damkohler  numbers,  on 
the  turbulent  flame  speed  in  the  form  (u  /q^)  and  on  the  orientation 
of  the  flame,  i.e. , on  B.  The  definitions  of  P and  K (cf.  Eqs.  (4.1) 
and  (4.  2)  ) yield 


(u  /q^).  = “ ^ (u  /q^  ) (P.  /P  ) 

o i o “ o i o 


tan  S X / P- 

o / I . I 


=-^  K 2 

‘ 1 + tan  0 p o 

o o 


,P.  K.  . 

. o 


(4.31) 


(4.32) 


for  i -•  1,  2 . 

Numerical  Results 

The  first-order  solutions  require  only  values  of  the  heat  release 

oa’-ameter  T and  of  c , and  m ; we  have  assumed  a range  of  values 
^ o 

for  T but  taken  C = 0.  3 and  m = 2 as  in  land  II.  For  one  case 

several  values  for  c are  assumed  to  demonstrate  the  relative 

o 

insensitivity  of  our  results  to  the  "cold  boundary"  condition.  The 
perturbation  solutions  on  the  other  hand  require  specification  of  other 

A 

parameters;  in  addition  to  the  eigenvalues  P^  and  the  angle  0^ 

appears  and  therefore  a value  for  ^ is  required.  Also  required  are 

values  for  c and  for  the  individual  components  of  m,  i.e.,  for  p 
mo 


Fiquation  (4.31)  is  based  on  the  assumption  that  the  factor 


7 

multiplying  (^  /u  ) in  Eq.  (4.  1)  is  invariant  with  respect  to  the 
0»  0 * 

perturbation  parameters  and  therefore  equal  to  ♦. 
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and  k.  Subsequent  applications  of  the  perturbation  solutions  with  these 

quantities  specified  require  values  for  C and  TTf  . Without 

1 o 

complete  rat'onali^.ation  we  take  p - k - 1 while  for  ♦ we  take  the 
strong  interaction  value  given  in  Appendix  3:  $ = 0.  1.  Generally  we  assunu 

calculation  the  effect  associated  with 
'^mo  ~ Note  that  in  Appendix  2 we  discuss  the  elements  of  laminar 

flame  theory  used  to  establish  these  values  which  are  close  to  the  value 

of  c = 0.  7 used  in  I and  II. 

mo 

The  asymptotic  approximations  to  the  first  order  solutions  in 

the  neighborhood  of  c = 1 are  applied  at  c = 0.95.  I(  is  found  necessary 

to  initiate  the  integration  of  the  perturbation  solutions  closer  to  c ~ 1 

in  order  to  achieve  satisfactory  accuracy;  accordingly,  we  take  ^ = 0.  98 

as  the  starting  po int  for  the  perturbation  solutions. 

For  purposes  of  exposition  we  consider  the  case  T = 5 , ^ = 0.  02 

’ o * ’ 

‘^mo  ~ to  be  prototypical  and  display  the  results  therefor.  In 

Fig.  2 we  show  the  first  order  solutions.  For  this  case  3 =87.2 

o ’ 

= 0,0897  and  the  angle  is  found  to  be  f).45°.  In  Figs.  3 and  4 

we  show  the  perturbation  solutions  for  Reynolds  and  Damkohler  numbers 

respectively.  For  the  former  we  take  Cj  = 1 , an  indicative  value. 

The  perturbation  eigenvalues  provide  the  results  of  greatest 

physical  significance  from  this  analysis.  In  Table  1 we  show  results 

for  several  cases,  i.e.,  for  several  values  of  t and  c^  and  c 

o mo  ‘ 

Presented  are  the  coefficients  a.^,b.^  (cf.  Eqs.  (4.  27)-(4.  30)  ) from 
which  the  ratios  (P./3^)  and  (K./K^)  may  be  computed.  Also  shown 
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are  the  ratios  of  the  physicallvr  slgaificant  quantities  (u  /q^)./(u  /q^) 

^ o ® i o ® o 

and  Eqs.  (4.  3 1)  and  (4.  32)  ) . For  the  f -perturbations  we  ha  ve 

again  taken  C ^ = 1 . 

Although  we  shall  discuss  the  implications  from  these  results  in 

detail  in  Section  6,  it  is  perhaps  appropriate  to  make  several  remarks. 

Note  that  apparently  the  greatest  sensitivity  of  the  turbulent  flame  speed 

and  of  the  orientation  of  the  flame  is  assoc  iated  with  the  influence  of 

molecular  transport,  i.  e.,  with  the  ^-perturbations.  The  extent  of 

this  influence  depends  on  the  heat  release,  decreasing  as  T decreases. 

In  fact  for  T = 5 it  appears  that  the  first  term  in  the  series  in  c 

-2 

cannot  be  continued  beyond  f = 10  without  raising  doubt  as  to  the 

-2 

validity  of  the  perturbation  analysis.  A value  of  c = 10  corresponds 
3 

to  j ^ standard  value  of  a,  namely  0.09,  is  assumed. 

Note  that  the  effect  of  finite  Damkohler  numbers  increases  as  r decreases 
and  that  the  perturbation  results  are  relatively  insensitive  to  the  value 
of  c 

mo 

The  influence  of  finite  Damkohler  numbers  on  flame  properties 

depends  directly  on  the  magnitude  of  the  parameter  However, 

4 

values  thereof  on  the  order  of  10  appear  to  be  required  in  order  for 
the  influence  of  finite  reaction  rate  to  be  comparable  to  that  of  finite 
Reynolds  numbers.  In  addition  the  K^-eigenvalues  are  exceedingly 

small.  Thus  the  changes  in  the  orientation  of  the  flame  due  to  finite 

• • ^ 

Damkohler  effects  are  due  essentially  to  the  0^ 


[ 

r 

5.  APPLICATION  TO  NORMAL  FLAMES 

We  now  apply  the  same  techniques  used  for  flames  undergoing 
strong  interaction  to  flames  involving  only  dilatation  due  to  heat  release, 
i.e.,  either  to  normal  flames,  or  to  unconflned,  oblique  flames. 
Because  the  treatment  closely  follows  that  in  the  previous  section,  we 
need  only  outline  the  steps  in  the  analysis.  In  this  case  the  turbulent 
kinetic  energy  is  conveniently  replaced  by  Q where  (1  - Q^)  Q t 

while  the  two  eigenvalues  to  be  determined  as  part  of  the  solution  are 
conveniently  taken  to  be 


(aCf,/(2c  - l)^,)(q  /u^)(l  - Q„) 

\ Im  2/  0 0 ® 

(5.1) 

A 

Q = 

Q /(I  - Q ) 

(5.2) 

The  quantities  (3  and  Q determine  the  turbulent  flame  speed  in  the 

form  (u  /q  ) and  the  ratio  of  the  turbulent  kinetic  energies  downstream 
o o 

and  upstream  of  the  flame,  i.e.,  . 

Equations  (3.  1)  - (3.4)  with  the  modelling  and  representation  of  the 
molecular  transport  terms  associated  with  Eqs.  (3. 6) -(3.  8)  provide 
the  starting  point  for  the  analysis.  With  9 = 90°  and  with  the 
alterations  of  the  equations  as  a result  of  the  different  eigenvalues, 

Eqs.  (4,  3)  - (4.  5)  are  changed  in  obvious  ways.  Finally,  if  the  new 
dependent  variables  given  by  Eqs.  (4.6)  are  introduced  and  again  c is 
treated  as  the  independent  variable,  we  find  the  following  first  order 


equations: 
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/i  + Q)c(l  -C 

V ^ CM  j.  rX 


— - ^ = (2c  - 1) 

dc  ^ ' dc  S(l+Tc)‘ 


£ ^ (1  +€C  + e (2  — (RE  S) 

c(l.c)'  ' ' \'l  ^ dc  2 


- 2(c_  - c)  ~ (RE,  S)  ) 
dc 


(5.3) 


dc 


D 

S 


dc 


G 

S 


A 

<T 

I + Tc 


(S  + Q)  + e 


d?  ' ? 


t^<t 


+ 2)S 


g)  (I 


+ TC  )' 


o 


(5.4) 


(5.5) 


dc 


G 

S 


(5.6) 


The  same  perturbation  parameters  C and  6 are  retained. 

These  equations  are  consistent  with  Eqs.  (4.  7) -(4.  10)  if  K = 0 and  if 
Q becomes  indefinitely  large:  note  that  Q~  -1  and  P if  Q ~ ® . 

flD  ^ ^ CD 


The  perturbation  equations 

The  expansions  given  by  Kqs.  (4.  11)  are  again  employed  with 

A 

those  for  (3  and  K replaced  by 

= K +,?,  t«3  .... 

0 12 

Q = Q +€Q,+6Q,+... 
o 1 2 
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When  all  of  the  expansions  are  substituted  into  Eqs.  (5.  3)  -(5.6), 
the  terms  independent  of  € and  6 yield 


dS  + Q )c(l  - c) 


0 = 1 °__ 

dc  S (1  + Tc  )' 

o 


(5.8) 


dG  G 

— ^ - ^ = C(3  +Q  ) — 

~ o o o , ~ 

dc  o 1 f Tc 


(5.9) 


dQ  G 
o o 


(5.  10) 


Equations  (5.  8)  - (5.  10)  are  essentially  those  treated  in  I with  the 
changes  due  to  redefined  eigenvalues  and  to  the  introduction  of  variable 
density  effects  via  the  exponent  m. 

Collection  of  the  first  order  terms  in  C and  6 from  Eqs.  (5.3)  - 
(5.6)  yield  the  equations  for  the  perturbation  solutions,  namely 


dD  D /uo.  . aa 

fi.,  ( ^-^)=-(2c  .1)(—  +11-  — 

iZ  \ S ^ mo  \ 

nr  n ' nn  rlr* 


dc  O 


Q.  +Q.  S. 

I — +— i ^ - 6 - 

S +Q  ^o  ^2  c 


o o o 


c (1  - c ) 


®ii'=i«o7-))  (2-r 

o'  dc 


o 2 o 


- 2(c  - c ) ~ (R  E.S  )) 

mo  ~ ' o 1 o / 


(5.11) 


dP 

o 

D 

o 

dc 

~ S 

o 

dG. 

1 

G 
o 

dc 

" S 

o 

VG 

o 

(5. 12) 


♦(1  + Tc  )T^(r  + 2)S 


o ' 


(5.  13) 


dQ.  G G S 

1 ^ _£  '^i  _i_  \ 

S G ~ S ' 

clc  o o o 


(5.  14) 


where  again  we  use  the  Kronecker  delta  for  compactness:  i = 1 for 
the  f-perturbations,  I = 2 for  the  6 -perturbations. 

The  considerations  made  previously  regarding  the  determination  of  the 
two  sets  of  perturbation  solutions  by  these  equations  apply  in  this  case. 
However,  in  place  of  Eq.  (4.  19)  we  find  for  the  function  P 


R„(Q„  + Q„)c(l  - c)  . 

-2—^ 2_t (♦M  (1-Q.  )) 

(l+Tc)*'  ' ° ^ 


(5.  15) 


For  convenience  we  introduce  M = ♦ M (1  - Q ) and  consider  it  to  be 

o ®o 

a parameter  in  the  perturbation  solutions. 


The  boundary  conditions 

The  previous  comments  concerning  specification  of  boundary 
conditions  at  c - c^  and  in  the  neighborhood  of  c 1 apply  here  as 
well.  We  should  note,  however,  that  as  shown  in  I the  results  for  first 
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order  solutions  for  normal  flames  depend  somewhat  on  the  value  of 

^ ; this  is  interpreted  as  indicating  the  dependence  of  turbulent  flame 
o 

speed  for  normal  flames  on  the  method  of  flame  stabilization.  Accordingly, 
we  expect  the  perturbation  eigenvalues  to  be  somewhat  dependent  on  c^  . 

The  solutions  for  the  range  0 ^ c c^  1 are  given  by  Eqs. 

(4.  20)  - (4.  22)  but  Eqs.  (4.  23)  - (4.  25)  are  replaced  by 

G = - (1  - Q ) + «T(1  + Q ) c (5.  16) 

0 o o 

G.  - Q.  = - G + CT  Q.  c,  i = 1 (5.  17) 

1 1 o i 

= €T  Q ^c',  i = 2 (5.  18) 

These  solutions  for  the  first  order  equations  applied  at  c = c^ 
provide  the  boundary  conditions  which  determine  the  first  order  eigen- 
values, F and  Q . Similarly,  those  applicable  to  the  perturbation 
o o 

, A 

equations  determine  the  perturbation  eigenvalues  p.  and  Q.,  i = 1,  2. 

The  boundary  conditions  for  the  solutions  in  the  neighborhood  of 
c = 1 are  developed  in  detail  in  Appendix  1 and  involve  no  arbitrariness 
if  the  solutions  are  to  have  acceptable  behavior  near  c = 1 . Thus 
integration  of  both  the  first  order  and  perturbation  solutions  is  initiated 
at  a value  of  c suitably  close  to  c = 1 . For  simplicity  we  repeat  our 
earlier  assumption  and  take  the  quantity  to  be  unity. 

The  form  of  the  perturbation  solutions  and  related  matters 

For  normal  flames  the  first  order  solutions  require  specification  of 
the  parameters  r,  f,  m and  c^  . Again  in  this  case  we  seek  the 


perturbation  solutions  in  a form  requiring  specification  of  a minimum 
number  of  additional  parameters  until  specific  values  for  the  perturbation 
eigenvalues  and  until  reconstituted  perturbation  solutions  are  desired. 
Thus,  for  example,  we  determine  the  solutions  in  the  form 

(^)  + (5.19) 


with  similar  forms  for  and  G.  (cf.  Eq.  (4.26)  and  comments 

applying  thereto).  When  these  forms  are  substituted  into  Eqs.  (4.21) 

and  (S.  17)  for  i 1 and  into  Eqs.  (4.22)  and  (S.  18)  for  i = 2 with 

, we  obtain  two  pairs  of  algebraic  equations  for  the  pairs  of 

eigenvalues  in  the  form  (K/?f  ) and  (Q./Q  ).  The  results  are  con- 

1 o 1 o 

veniently  arranged  as  in  Eqs.  (4.  27)  - (4.  30) , namely  as 


(5.20) 


(5.21) 


V^O=“'>21 


(5.22) 


°2^0o=“'’22 


(5.23) 


where  the  a_  and  b_  coefficients  depend  only  on  the  parameters  determining 
the  first  order  solutions  and  indicate  the  sensitivity  of  the  perturbation 
eigenvalues  to  molecular  transport  and  to  the  parameters  Cj  and  M. 

The  physically  interesting  results  from  this  analysis  pertain  to  the 
perturbed  values  for  the  turbulent  flame  speed  and  for  the  ratio  of  turbulent 
kinetic  energy  downstream  and  upstream  of  the  flame.  These  are  found 


37 


from  Eqs.  (5.  1),  (5.  2)  and  (5.  7)  with  ^ taken  as  a constant  to  be 
expressible  in  terms  of  the  perturbation  eigenvalues  as 


(5.24) 


(5.25) 


Numerical  Results 

We  have  obtained  numerical  solutions  for  normal  flames  based  on 

our  standard  values  of  € = 0.  3 , m = 2 for  a range  of  values  of  the  heat 

release  parameter  T.  Generally,  we  take  c^  - 0.02  but  for  one  case 

we  show  the  influence  of  different  values  thereof.  As  in  the  strong 

interaction  case  the  perturbation  solutions  require  specification  of 

additional  parameters:  generally,  we  assume  c - 0.  833  but  again 

for  one  case  we  show  the  effect  of  c 0.75.  The  value  of  ^ is  assvimed 

mo 

to  be  1.4  (cf.  Appendix  3).  Finally,  we  again  take  p = k = 1. 

The  numerical  solution  of  the  equations  for  normal  flames  is 
found  to  be  increasingly  difficult  as  T increases.  With  the  standard 
method  for  solving  differential  equations  we  use, the  permitted  step-size 
becomes  exceedingly  small  as  T increases.  Thus  practical  considerations 
limit  our  results  to  T s 4 and  to  the  requirement  that  the  asymptotic 
solutions  in  the  neighborhood  of  c^  ~ 1 be  applied  at  ^ = 0.  95  for 
both  the  first-order  and  perturbation  solutions. 
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"We  take  the  case  of  r - 4 , c =0.02  and  c = 0.833  to  be 

o mo 

prototypical  and  show  the  results  therefor  in  detail.  In  Fig.  5 the  first- 
order  solutions  are  given  while  Figs.  6 and  7 give  the  perturbation 
solutions  associated  with  finite  Reynolds  and  Damkohler  numbers 
respectively.  Again  for  the  former  we  use  the  representative  value, 

Cj  = 1.  Note  from  Figs.  6 and  7 the  dominance  of  the  perturbed 
turbulent  kinetic  energy,  i.e.,  of  and  , over  the  other 

perturbation  functions. 

In  Table  2 we  present  for  a range  of  solutions  the  results  for  the 

eigenvalues  and  for  parameters  of  interest  based  thereon.  The  first 

order  eigenvalues,  the  coeffic ients  a.,  and  b..  which  permit  the 

U ij 

perturbed  eigenvalues  to  be  computed,  and  the  physically  significant 

quantities  based  on  the  perturbed  eigenvalues  are  given.  Again  we 

shall  discuss  these  results  in  detail  in  Section  6 but  certain  general 

remarks  can  be  made  here.  As  in  the  strong  interaction  case  the 

magnitude  of  the  perturbations  associated  with  large  but  finite  Reynolds 

numbers  increase  with  heat  release  but  decrease  with  corresponding 

values  for  Damkohler  numbers.  Moreover,  the  effect  of  the  perturbations 

on  turbulent  flame  speed  are  of  opposite  sign  but  on  the  ratio  of 

turbulent  kinetic  energies  are  of  the  same  sign.  Finally,  we  note 

the  weak  dependence  of  the  predicted  behavior  with  c and  c 

o mo  * 
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6.  INTRODUC9  ION  OF  LAMINAR  FLAME  PARAMETERS  AND 
COMPARISON  WITH  EXPERIMENT 

As  indicated  earlier  the  behavior  of  turbulent  flames  is  frequently 
correlated  by  exper irnontalists  in  terms  of  the  characteristics  of  the 
classical  laminar  flame,  determined  either  experimentally  or  theoretically, 
with  the  same  chemical  and  thermodynamic  situation  as  prevails  for 
the  turbulent  flame  in  question.  The  implication  from  this  practice 
is  that  one  or  both  of  our  perturbation  parameters  are  expressible  in 
terms  of  laminar  flame  characteristics. 

To  pursue  this  implication  further  we  consider  the  laminar  flame 
speed,  u^  , and  a measure  of  the  laminar  flame  thickness,  : from 
[11,  12]  we  have 


u,  K,(U  w Ipy 

I 1 o max  o 


A = K (1/  p /w  ) 
L Zoo  max 


1) 


where  Kj  and  are  non-dimensional  quantities  given  by  the 
appropriate  laminar  flame  solution.  In  the  present  context  they  should 
be  treated  as  functions  of  the  heat  release  parameter  T and  as  dependent 
on  the  chemical  kinetic  model  representing  w/w  . 

rxiaLX 

The  quantities  u^  and  given  by  Eqs.  (6.  1)  and  (6.2)  can  be 

readily  expressed  in  terms  of  our  perturbation  parameters:  we  find 
from  the  definitions  of  € and  6 


u^  = q (K^/aL)(€/6) 
jt  o 1 3 


(6.3) 


r 
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= ji^(K^a^I  ) €6 
X o Z 3 


(6.4) 


Equation  (6.3)  indicates  that  is  undefined  in  the  first-order  theory 

since  f,  6 = 0 while  Eq.  (6.4)  is  consistent  with  the  basic  notion  of 
the  Bray-Moss  model  regarding  a thin  laminar  flame  provided  both 
the  turbulent  Reynolds  and  Damkohler  numbers  are  large. 

Alternatively,  Eqs.  (6.3)  and  (6.4)  yield 


.3  ^ 

'''I  ''  ^ 


4 


KK^a 


K, 


(6.5) 


L_ 

/ u 3/2 

o i K^La 
2 3 


(6.6) 


which  show  that  our  perturbation  parameters  are  expressible  in  terms 
of  the  laminar  flame  characteristics. 

Frequently  only  one  of  the  two  laminar  flame  characteristics 
introduced  here  are  used  by  experimentalists  to  correlate  their  results, 
in  this  case  only  one  of  our  perturbation  parameters  can  be  eliminated. 
We  shall  illustrate  this  in  detail  later. 

Strong  Interaction 

In  the  case  of  strong  interaction  the  orientation  of  the  reaction 
/.one  is  the  predicted  quantity  most  directly  comparable  with  experiment. 
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From  Eqs.  (4.  27)-(4.  30)  and  (4.  32)  we  find 


tan  9 / 

1 + tan  9 \ 


or  for  0 ^<  \ 

o 


= ( 1 - i f ( ) f + ( ) 57!  e ' + . . . ) 

o y V O'  I 


(6.7) 


The  results  given  in  Table  1 can  be  used  to  evaluate  the  coefficients 
in  Eq.  (6.  7). 

The  value  of  «!'  has  been  selected  as  discussed  in  Appendix  3 so 
that  the  predictions  of  the  first-order  theory  for  0^  are  in  good  agree- 
ment with  the  results  of  Wright  and  Zukoski  [2],  There  are,  however, 
no  data  permitting  assessment  of  the  perturbation  effects.  In  this 
situation  it  may  be  useful  to  consider  the  counterpart  for  strong  interaction 
of  the  correlation  of  the  properties  of  normal  flames  given  in  [1]. 
Accordingly,  we  introduce  the  laminar  flame  speed  and  eliminate  the 

Damkohler  parameter;  from  Eq.  (6.  1)  and  the  definition  of  ^ we  have 

o 


V ’*'7r)  -i  *^1' 

so  that  Eq.  (6.  7)  becomes 


♦ M K,  q ^ 

o 1 ^»  \ 


^ ^*^21 ‘‘‘^22^  I 2 / 

3 


+ . • • 


(6.8) 


We  eliminate  the  turbulent  kinetic  energy  downstream  of  the 

flame  from  Eq.  (6.8),  since  this  is  generally  not  measured.  For  this 

purpose  only  the  first-order  estimate  of  is  needed;  from  the  ratio 

u /q  , the  orientation  of  the  flame  and  the  resultant  velocity  upstream 
o ^ ® 

of  the  flame,  V , we  have  for  6 1 

o o 


-2 

2 


P e V 

o o o 


(6.9) 


Thus  Eq.  (6.8)  becomes 


*">21  *'>22' 


M K p ^ 

Q 1 o o o', 


(6.  10) 
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In  Appendix  Z we  employ  laminar  flame  theory  to  provide  inter 

* 

alia  estimates  for  Kj  and  for  the  quotient  . To  illustrate 

the  results  given  by  Eq.  (6.  8)  we  take  the  exponent  n = 5 , C j =1 

and  T = 5 and  the  results  for  our  prototypical  case  T = 5,  c = 0.02. 

o 

With  these  values  Eq.  (6.8)  yields 

e~e  ( 1 + € fll.7  - 1.02(10"^)(V^/u^)  +...)  (6.11) 

o y \ o i / I 

From  Eq.  (6.  11)  we  see  that  for  values  of  the  velocity  ratio  V /u.  of 

o £ 

30-40,  corresponding  to  values  of  of  3-4,  the  term  arising 

from  the  replacement  of  the  Damkohler  perturbation  is  comparable  to 
that  due  to  the  finite  Reynolds  number  perturbation  but  of  opposite  sign 
so  that  the  insensitivity  of  the  angle  of  the  flame  observed  by  Wright 
and  Zukoski  [2]  to  a wide  variety  of  conditions  is  perhaps  explicable. 

For  a fixed  value  of  the  ratio  V /u,  less  than  34  a decrease  in  the 
turbulence  Reynolds  number  leads  to  an  increase  in  the  angle  of 
orientation;  the  opposite  is  true  for  values  of  the  same  velocity  ratio 

It  should  be  noted  that  two  separate  calculations  of  these  quantities 
can  be  envisaged.  The  origin  of  Kj  via  Eq.  (6.  1)  implies  that  it  is 
determined  by  the  solution  for  the  classical  laminar  flame;  on  the  contrary, 
the  quotient  M^/I^  relates  to  the  distribution  f(c)  within  a laminar  flame 
embedded  in  the  turbulent  reaction  zone  in  question  and  thus  in  principle 
involves  a separate  calculation.  For  simplicity  and  consistency  with  the 
assumed  insensitivity  in  the  Bray-Moss  model  of  the  results  to  the  details 
of  f(c),  we  use  the  same  solution  for  both  calculations. 
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1 

I 

greater  than  34.  Note  that  in  general  the  validity  of  the  perturbation 

-2  3 ' 

analysis  may  be  questioned  for  cJ®  10  , i.  e.,  for  < 10  , | 

Normal  Flames 

In  considering  the  implications  of  our  perturbation  analysis  for 
normal  flames  we  focus  on  the  recent  correlation  of  [ l]  which 
relates  the  ratio  of  turbulent  to  laminar  flame  speeds  to  the  ratio  of 
initial  turbulent  intensity  to  laminar  flame  speed  for  a wide  range  of 
experimental  data  and  thus  of  turbulence  Reynolds  numbers.  The  absence 
in  the  correlation  of  heat  release  as  a parameter  implies  that  over  the 
range  of  T of  practical  interest,  i.e.,  4 T < 9 , the  variation  of  the 

ratios  in  question  is  within  the  scatter  of  the  experimental  results. 

1 

Accordingly,  we  treat  our  results  for  T = 4 as  being  representative  of  j 
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Equation  (6.  12)  leads  to  the  anomalous  result  that  as  ^^/u  0, 

o ^ 

u /u  ~ 0;  this  pathological  behavior  is  associated  with  the  assiimplion 

of  a high  turbulence  Reynolds  number  inherent  in  the  Ilray-Moss  model 

and  can  be  removed  by  heuristically  replacing  the  left  side  of  Eq.  (6.  12) 

and  appropriate  subsequent  equations  by  (u  /u  ) -1  . 

O X 

Equation  (6,  12)  with  ( = 6=0  and  with  the  left  side  altered  provides 
the  basis  on  which  the  value  of  ^ = 1,4  used  here  is  selected  in  Appendix  3, 
To  investigate  the  effect  of  Reynolds  number  given  by  Eq.  (6.  12)  we 

A 

eliminate  M and  introduce  the  laminar  flame  speed.  From  the  definition 

A 

of  M and  6 we  find  that 

q. 


M6  = 


'o  N 


A \ 7 ^ 

3 1 + Q n' 

o I 


(6.  13) 


so  that  Eq.  (6.  12)  becomes 


u 

- 1 = ( 


|5  (1  +Q  ) 
o o 


4r  a ^ 
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^ ^ V^o.,o^^l2'*^^l  ^12^  ■‘■^^11  '*■^1^11' 


^ <Q-o"22-^‘^2 


2.) — \ — I2. 

^ ^ ^3  ^ 1 + Q u?  f 


o I 


(6.  14) 


To  evaluate  the  coefficients  in  Eq.  (6.  14)  we  use  the  results  in 
Appendix  2 and  of  our  protypical  case  for  normal  flames;  thus  we  let: 

T = 4,  c^  = 0.02,  ♦ = 1.4,  Cj  = 1 , and  n = 5,  There  results 

Uq  q 

^ - 1 = 1.14  ^ (1  +(11.3  . 0.750  -£)^+...)  (6.15 

/ £ u^ 
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Equation  (6.  15)  displays  the  same  opposing  effects  as  shown  in  Eq.  (6.  11) 

for  the  strong  interaction  case;  thus,  for  q^/u.  less  than  3.  88  a decrease 

o X 

in  turbulence  Reynolds  number  increases  the  turbulent  flame  speed 

whereas  for  q^  /u^  greater  than  3.  88  the  opposite  effect  is  predicted. 

Figure  8 shows  u^/u^  calculated  from  Eq.  (6.  15)  for  € = 0, 

-3  -3 

3 (10  ) and  7 (10  ).  In  the  region  q*/u,  < 3.  88,  the  increase  in  the 

O A 

ordinate  resulting  from  these  values  of  C is  too  small  to  show  in  this 
figure.  Consequently,  the  perturbation  analysis  predicts  that,  at  a 
fixed  turbulence  Reynolds  number,  the  ratio  of  turbulent  to  laminar 
flame  speeds  initially  increases  nearly  linearly  with  increase  in  upstream 
turbulence  intensity.  At  higher  turbulence  intensities  the  finite  Reynolds 
number  curves  drop  below  this  straight  line.  The  curves  for  € 0 

in  Fig.  8 are  continued  until  the  perturbation  reaches  20%. 

We  note  that  the  predictions  of  the  theory  similar  to  Eq.  (6.  15)  but 
for  different  values  of  t and  of  the  exponent  n in  the  laminar  flame 
theory  do  not  appear  to  differ  significantly  in  a qualitative  sense  from 
those  given  here. 

Figure  9 shows  two  sets  of  experimental  data  which  approximately 
define  two  extremes  in  the  correlation  of  Abdel-Gayed  and  Bradley  [l]. 

The  upper  band  is  from  [14]  and  the  lower  band  from  [15],  both  being 
interpreted  in  the  manner  of  [l].  It  has  been  assumed  that  the  length 
scale  in  our  turbulence  Reynolds  number  R^  is  the  same  as  that 
used  in  [ l]  to  define  R^  = u'  l/v , where  u'  is  the  root  mean  square 
of  the  x-wise  velocity  component.  In  converting  from  the  abscissa 
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2 — ~ 

u /n'  of  [ l]  it  has  been  assumed  that  pu*  /pq  = 0,3  and  that 
JL 

i -3  -3 

R = R./O.  3®.  Thus  f = 0,  3 (10  ) and  7 (10  ) corresponds  to 

T I 

3 2 

R 00  and  R = 2 X 10  , 8.  8 X 10  , respectively. 

JL  JL 

If  the  two  sets  of  data  are  given  equal  weight  then,  superficially 
at  least,  a trend  suggesting  a reduction  of  turbulent  flame  speed  at 
lower  turbulence  Reynolds  numbers  might  be  discerned,  as  proposed 
in[l].  Closer  inspection  of  the  data  in  Figure  9 argues  however  that 
no  convincing  Reynolds  number  trend  is  evident  in  the  data  of  either 
[ 14]  or  [is],  considered  separately,  despite  the  fact  that  each  covers 
a significant  Reynolds  number  range.  There  is  clearly  serious  dis- 
agreement between  data  from  these  sources  at  comparable  Reynolds 
numbers  and  the  principal  distinction  would  appear  to  relate  to 
experimental  configuration  and  technique.  The  emphasis  of  these 
remarks  is  not  significantly  changed  if  the  other  data  from  [ l]  is 
incorporated  in  Figure  9. 

A comparison  between  the  perturbation  solutions  and  experiment 
is  therefore  not  straightforward.  If  experimental  data  are  chosen  from 
a particular  source  such  as  [l4],  then  the  predicted  trend  with  Reynolds 
number  is  not  supported.  The  relatively  narrow  range  of  flame  speeds 
within  which  the  perturbation  is  valid  is  comparable  to  the  scatter  in 
the  data  [14],  On  the  other  hand  the  trend  predicted  by  the  perturbation 
analysis  agrees  with  that  reported  in  [l].  It  is  important  to  note  that 
the  form  of  presentation  of  the  experimental  data  adopted  here  is  that 
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suggested  most  naturally  by  the  present  analysis  and  not  that  of 
Abdel-Gayed  and  Bradley  [l].  Their  abscissa,  u /u^  , is  related 
to  the  inverse  of  that  introduced  in  Fig.  9 and  leads  to  curves  of 
broadly  hyperbolic  character.  In  consequence,  and  in  contrast  to 


the  present  approach,  attention  is  focussed  on  modest  values  of 


The  Influence  of  the  Perturbations  on  the  Structure  of  the  Normal  Flame 
It  is  of  interest  to  consider  the  influence  of  the  perturbations 


on  (he  structure  of  the  turl)ulent  reaction  zone.  To  show  this  influence 
requires  selection  of  a length  scale  model.  For  consistency  we  (akc 
q i to  be  constant.  In  this  case  the  x-coordinate  is  given  by 
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We  see  from  Eq.  (6.  16)  that  the  predicted  spatial  distributions  within 

the  reaction  zone  depend  on  both  the  perturbed  eigenvalues  via  the  ratios 

(u  /q  ).,  i = 1,  2 and  on  the  perturbed  solutions  themselves  via  S.(c). 
o o 1 1 

In  the  spirit  of  our  earlier  considerations  we  eliminate  the  Damkohler 
perturbation  in  favor  of  the  laminar  flame  speed;  there  results 


(u  /q") 
o o o 


1 

o ^ o 1 


(u  /q  ) 
o o o 


, ,M  5 ,j,  (U  /q=), 

11/  O ^ ^ Xa 

3 u‘^  M ) ^r. 

6 o 0^000 


J + 


(6.  17) 


If  we  evaluate  the  coefficients  in  Eq.  (6.  17)  on  the  basis  of  onr 
prototypical  case  for  normal  flames,  we  find 


- I 

— = 0,079ll 

o \ 


J - e;  J,  + 11.2 J + 0.  556  — 
o V 1 o u. 

I 


( - 2.  55  J 

M 

o 


o)  I + . . . ) 


(6. 18) 


We  again  see  the  competing  effects  of  the  Reynolds  and  Damkohler 
perturbations.  In  Fig. 10  we  show  the  distribution  of  the  mean  product 
concentration  according  to  the  first-order  theory  and  to  Eq.  (6.  18)  with 
q^/u  = 10  which  corresponds  to  the  limit  of  validity  of  the  perturbation 
analysis  for  the  assumed  value  of  the  perturbation  parameter, 

C=  3(10‘^)  . 
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Several  remarks  regarding  Fig.  10  are  indicated.  First  note 

that  the  thickness  of  the  reaction  zone  is  on  the  order  of  the  reference 

length  scale  I ; this  is  in  accord  with  the  results  of  I but  we  here 
o 

explicitly  incorporate  density  effects  in  the  modelling  and  therefore 
have  somewhat  altered  solutions.  The  effect  of  the  combined  Reynolds 
and  Damkbhler  perturbations  on  the  structure  of  the  reaction  zone 
tends  to  reduce  the  thickness  for  the  particular  value  of  the  ratio 
q selected.  For  smaller  values  of  this  ratio  the  predicted  reaction 

zone  would  be  thickened  somewhat. 

7.  CONCLUDING  REMARKS 

Premixed  turbulent  combustion  theory  [5-8]  has  been  extended  to 
include  effects  of  molecular  transport  and  finite  chemical  reaction  rates. 
The  corresponding  premixed  laminar  flame  properties  then  arise  quite 
naturally  as  parameters  of  the  turbulent  flame  formulation,  as  already 
found  in  experiments.  It  is  the  first  time  that  this  has  been  demonstrated 
theoretically. 

A perturbation  analysis  has  been  made,  for  large  values  of 
turbulence  Reynolds  number  and  Damkohler  number,  and  applied  to 
highly  oblique  confined  flames,  in  the  strong  interaction  limit,  and 
to  normal  or  unconfined  oblique  flames.  Numerical  solutions  have 
been  obtained  using  separate  values  of  the  important  modelling  parameter 
4 for  these  two  cases. 
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The  strong  interaction  solution  yields  an  expression  (Eq.  (6.  11)  ) 
for  the  flame  angle  6 as  a function  of  and  I*  i®  observed 

that,  in  the  range  of  validity  of  the  perturbation,  the  perturbation  terms 
arising  from  Da  and  are  of  similar  magnitude  but  opposite  sign. 

The  tendency  for  these  terms  to  cancel  each  other  is  offered  as  a 
tentative  explanation  of  the  experimentally  observed  [2]  insensitivity 
of  the  flame  angle  to  large  changes  in  the  composition  and  temperature 
of  the  combustible  mixture. 

In  the  case  of  normal  and  unconfined  oblique  flames  the  turbulent 

flame  speed  in  the  form  u ®/u  . has  been  predicted  as  a function  of 

o * 

turbulence  intensity,  q /u  , and  turbulence  Reynolds  number,  R~  ; 

O JL  ^ 

see  Eq.  (6.  15).  These  are  essentially  the  same  variables  as  those 
employed  in  an  empirical  correlation  of  experimental  data  by  Abdel- 
Gayed  and  Bradley  [l].  The  perturbation  solution  predicts  that  the 
turbulent  flame  speed  at  finite  R^  falls  below  the  high  Reynolds 
number  limit  by  an  amount  which  increases  with  increasing 
and  with  decreasing  • The  analysis  shows  that  this  effect  arises 
because  of  the  finite  reaction  rate  terms  in  the  equation,  rather  than 
the  molecular  transport  terms,  whose  direct  influence  on  u^/u^  is 
negligible.  Equation  (6.  15)  then  shows  that  the  first-order  approximation 
to  the  flame  speed  is  valid  so  long  as 
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Convincing  support  for  the  perturbation  predictions  for  normal 
flames  is  not  found  in  the  experimental  data.  Individual  experiments 
appear  not  to  show  a trend  of  variation  with  although  the 

scatter  of  data  points  is  large.  On  the  other  hand  the  predicted  trend 
agrees  qualitatively  with  that  found  by  Abdel-Gayed  and  Bradley  [ l] 
from  their  correlation  of  experimental  data. 
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?*/«,  = 10 
O I 


APPENDIX  1.  I ho  llohii  vior  as  > - I 

Wc  consider  the  behavior  of  the  first  order  and  perturbations 
solutions  in  the  neighborhood  of  c = 1 , first  for  the  case  of  strong  inter- 
action and  subsequently  for  normal  flames.  If  Eq.  (4.  12)  is  put  in  an 
approximate  form  appropriate  for  c~l,  it  may  be  solved  to  yield 


S = -X(l  - c) 
o 


(Al.  1) 


where  X 


= *('-( 


1 + 4P  /(I  + T) 
o 


Next  the  right  side  of 


Eq.  (4.  13)  can  be  similarly  approximated  so  that  with  Eq.  (Al.  1) 


integration  leads  to 


G = - aX(l  - c)/(l  - X) 
o 


(A1.2) 


rhere  a=|^<T/(l  + T)^  Eqs.  (A  1 . 1)  and  (A  1. 2)  permit 


Eq.  (4.  14)  to  be  integrated  so  that 


Q = -a(l  - c)/(l  - X) 
o 


(A1.3) 


Equations  (A1)-(A3)  applied  at  a value  of  c in  the  neighborhood  of 
c = 1 provide  the  boundary  conditions  for  the  first  order  solution  and 
are  employed  in  I and  II.  Note  that  these  solutions  are  free  of  arbitrary 
constants,  the  parameters  X and  a depending  only  on  the  eigenvalues 

A 

(3^  and  . In  the  solution  for  a second  root  of  the  secular  equation 


for  X is  discarded  so  that  ~ 0 as  c ~ 1 ; furthermore,  in  the 


integrations  leading  to  Eqs.  (A  1.2)  and  (A  1.3)  complementary  solutions 


which  have  unacceptable  behavior  as  c ~ 1 must  be  suppresaed, 
leaving  the  particular  solutions  as  the  only  valid  ones.  Similar 
considerations  arise  with  the  perturbation  solutions. 


As  a preliminary  remark  we  recall  that  earlier  we  assumed  the 
parameter  to  be  unity;  in  addition  we  note  that  the  quantities 
Ej(^)  and  ^^(c)  and  the  ratio  (V Iv^  will  approach  constant  values 
as  c ~ 1 : we  shall  denote  these  values  with  the  subscript  “.  Finally, 
the  asymptotic  solutions  for  the  first  order  functions  given  by  Eqs. 

(A  1.  1)  - (A  1. 3)  are  used  in  the  analysis  of  the  asymptotic  behavior  of  the 
perturbation  solutions. 

Equation  (4.  15)  takes  on  the  following  form  as  "c  ~ 1 : 


dS.  S. 

— - + — — = - (1  - X)  K. 

dc  1 . c 


(A1.4) 


where  K. 


P-  2(c  ^ - 1)  , V 

K.  = - 6.,  , i i»l  + 6.,C  — 

‘ p ‘2  - I „ o >11  K, 

o 


N / 2E, 

X / 2« 


2(c  - 1) 

mo 


il  1 - X \ 2c  - 1 2c 

mo  mo 


— E,.) 


We  note  that  again  the  complementary  solution  must  be  suppressed; 
thus  the  solution  to  Eq.  (A  1.4)  is 


■ ^ 

j 


(A1.5) 


S.  = - 

I 


X(1  - A.)  K.  (1  - c) 


1 - 2X 


With  determined,  the  asymptotic  approximation  to  Eq.  (4.  17)  yields 


dG.  , G. 
‘ +-!■  ‘ 


aK 


X ^ 
dc  1 - c 


i , 2 / i \ - aX  , « V 2 +p 

1 - 2X  '^o’^  \lC  / ■ ®il  I - X ^ 


= L. 


(A1.6) 


with  the  solution 


XL.(1  - c) 

^’i  " TTx 


(A1.7) 


Finally,  a convenient  form  for  the  boundary  conditions  for  the 
functions  is  obtained  by  substitution  of  Eq.  (4.  18)  into  Eq.  (4.  17), 
by  applying  the  approximations  appropriate  for  1 , and  by 
quadrature;  the  result  is 


v"(ir)<'  - S'! 


{A1.8) 


Equations  (A  1 . 5),  (A  1.7),  and  (A  1 . 8)  applied  at  a value  of  c in  the 
neighborhood  of  c = 1 provide  the  appropriate  boundary  conditions  to 


initiate  an  integration  in  the  direction  of  ^ for  the  strong  interaction 


case. 


1 


Equations  (Al,  1)-(A1*3)  apply  for  the  first  order  equations  for 
normal  flames  provided 


X = i ^1  - (l  + 415^6  ^/(l  j 


a = - € T Ci  /(I  +r) 

O 


(AI.9) 


(Al.  10) 


In  addition  the  solutions  for  the  perturbation  functions  S and  G 

i i 

given  by  Eqs.  (A  1 . 5)  and  (A  1 . 7)  are  retained  with 


Q. 

-i-  + 

^ Q 


2(c  - 1)  Q V 

„ i . I , mo  o ® 

K,  = — + — -6..  ; rM+a,C,  — 

i2  2c  - 1 k il  1 1/ 

mo  (1  + T)  o 

o o 


. , 2E  2(c  - 1) 

6 ^ i ^ mo 

ill-X>2c  -1  2c  “1  1®/' 

mo  mo 


(Al.  11) 


orK  €tQ  ,Q.  , , 

o 


(Al.  12) 


Finally,  a relation  for  the  perturbations  convenient  for  numerical 
analysis  and  completely  analogous  to  Eq.  (A  1.8)  is  obtained  from  Eqs. 

(5.  13)  and  (5.  14)  with  approximations  appropriate  for  1,  i.e., 

' ^i  = ^i  ^®il<^  ^’■>^^*’^0  (A1.13) 

^o 

W»!  nt>te  that  these  results  for  the  asymptotic  behavior  of  the 
perturbation  solutions  as  c ~ 1 for  normal  flames  are  consistent  with 


those  for  the  strong  interaction  case  if^Q~^,  Q 1, 

o o o o 

A A Q 

K — 0 , Q M ~ 771  > replacements  which  correspond  to  0 ^ 90  and 
o o 


Q becoming  indefinitely  large. 


APPENDIX  Z.  Determination  of  Perturbation  Parameters  by 


Laminar  Flame  Theory 

In  earlier  studies,  Bray  and  Moss  [5,6]  evaluate  the  integrals 

^k  + 3 " S k = 0,  1 (A2.  1) 

0 max 


which  arise  in  the  closure  of  the  chemical  source  term,  on  an  ad  hoc 
basis.  The  burning  mode  pdf,  f(c),  is  taken  to  be  a simple  battlement 
shape,  broadly  characteristic  of  a laminar  premixed  flame,  while 
the  reaction  rate  employs  an  Arrhenius  form.  The  requirements  of 
the  perturbation  analysis  for  similar  information  are  greater  and 
invite  a more  rational  approach. 

The  species  balance  equation  for  c in  a one-dimensional  laminar 
flame  may  be  written  in  the  form 


Pu  = m = constant  (A2.3) 

If  we  assume  the  Schmidt  Number  to  be  unity  and  introduce  the 
simplifying  thermodynamic  description  of  Section  2,  then 


HI  = PD 


I 

t 


■ 


whence  Eq.  (A2.  2)  becomes 


(14.  — 1 = 

dx  dx  1.  ^ ' dx  J 


(A2.5) 


If  we  introduce 


S = u 


(1  + Tc)  ^ 
m dx 


(A2.6) 


and  treat  c as  the  independent  variable,  then  Eq.  (A2.6)  may  be 


written  as 


d<; 

S(l  - ~)  = — - (1  + Tc)  w(c) 
m 


Following  Bray  and  Moss  [5,6],  we  have 


(AZ.7) 


= it  I 73^ 


dx  Jq  (dc/dx) 


(A2.8) 


or,  from  Eq.  (A2.6), 


f(c)  = 


1 + T C 

S \ (d  + Tc)/S  ) dc 
0 


(A2.9) 


We  note  that  the  integrals  indicated  in  Eqs.  (A2.8)  and  (A2.9)  have  been 

implicitly  taken  to  be  convergent.  In  general  this  is  not  the  case  and  a 

concept  analogous  to  that  used  in  experimental  turbulence,  namely  gating, 

must  be  introduced.  The  consequence  is  that  f(c)  is  defined  only  for 

the  range  C ^ c * 1 - C where  0 < c 1 : entries  for  0 s c ^ C 
g g g g 


r 

t 

I 

and  for  1 - c 1 are  attributed  to  the  Dirac  delta  functions  at 
c = 0,  1.  However,  we  shall  see  that  when  f(c)  is  used  in  evaluating 
1 the  various  integrals  of  interest,  the  integrals  are  convergent,  the 

integral  in  the  denominator  of  Eq.  (A2.  9)  is  inessential,  and  c 

g 

can  be  allowed  to  vanish.  Accordingly,  the  determination  of  S(c) 
for  a specified  reaction  rate  expression,  w(c),  permits  an  entirely 
self-consistent  determination  of  the  moments  and  of  other  quantities 
required  in  the  perturbation  analysis. 

Spalding  [14]  identifies  certain  reaction  rate  expressions  which  are 
both  plausible  and  permit  analytical  solution  of  Eq.  (A2.7)  subject  to 
the  appropriate  boundary  conditions.  It  is  thus  convenient  and  sufficient 


S(0)  = S{1)  = 0 . 


(A2. 12) 


By  inspection  the  solution  is  found  to  be  especially  simple,  namely 


I 


S = c(l  - c"‘S 


X = n 


(A2.  13) 


so  that  Eq.  (A2.9)  becomes 


f(c)  = 


1 + T c 


-1 


c(l  - c ) ^ ( 1 + Tc)  ^c  ( 1 - c"^  S dc 


€ c s 1 - f 

g g 


(A2.  14) 


where  we  explicitly  introduce  the  gating  parameter,  € 


We  now  use  these  solutions  to  obtain  with  € ~ 0 the  parameters 

g 


'mo  1 

^ c dc 
0 


n + 1 


(A2.  15) 


o _ H \ ' J - c ■, 

S ® *^0  1 - c"'^ 


(A2.  16) 


Note  that  c is  independent  of  j but  that  (M  /!,)  does  depend  on 
mo  0 3*^ 

T via  O'.  In  general  a numerical  evaluation  of  Eq.  (A2.  16)  is 


required.  In  addition  the  parameters  and  arising  in  Section  6 


can  be  evaluated.  First,  by  comparison  of  the  definitions  of 
and  X (cf.  Eqs.  (6.  1)  and  (A2.  11))  we  find 


APPENDIX  3.  Numerical  Values  of  O 


The  parameter 


aC  i 


4> 


1 


(2c  - 1)  I 

mo  o 


(A3.  1) 


must  be  specified  before  the  perturbation  solutions  can  be  calculated 
and  before  comparison  with  experiment  can  be  carried  out.  In  the 
strong  interaction  case,  ♦ is  estimated  from  the  experiments  of 
Wright  and  Zukoski  [2]  by  matching  the  measured  flame  angle  to  that 
predicted  by  the  first-order  solution.  The  relevant  equation  is  (see 
Eq.  (4.1),  (4.2)) 


2 

tan  6 


(A3.  2) 


0 (t 
o o 

If  we  use  our  prototypical  case  for  strong  interaction  to  select  a value 

A 

for  $,  we  find  with  B = 87.2,  K = 0.0897  and  $ = 0.  1 that 

o o 

= 6.45°  in  satisfactory  agreement  with  [2].  Accordingly,  we 
adopted  this  value  of  ♦ for  all  calculations  of  the  limiting  case  of 
strong  interaction. 

For  normal  flames  and  unconfined  oblique  flames,  4 is 
estimated  from  the  experimental  data  of  [141.  Using  the  first-order 
solution,  we  obtain  $ from  Eq.  (6.  12)  as 


u 


^ - 1 = r — ? — r 

(1  + Q ) -^  "i 
^o  o 


(A3.  3) 


With  B = 0.413  and  Q = 1.62,  it  is  found  that  $ = 1.4.  In  Fig.  9, 

the  line  f = 0,  representing  the  first-order  solution  with  = 1.4, 

is  seen  to  compare  well  with  experimental  data  from  [14]. 

It  is  evident  from  (A3.  1)  that  the  parameter  $ incorporates 

uncertain  features  of  both  the  thermochemical  and  fluid  mechanical 

modelling.  The  significantly  different  values  of  $ necessary  to 

admit  quantitative  comparison  in  both  the  unconfined  and  strong 

interaction  cases  suggest  major  differences,  either  in  local  structure, 

f(c)  (and  hence  c ),  or  in  the  "universal"  parameters,  for  example, 
m 

On  the  other  hand,  in  the  light  of  earlier  discussions,  the  experimental 
data  are  also  not  without  ambiguity.  Further  investigation  is  indicated. 
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